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Abstract 



The object of this paper is to present a computable solution of a fractional partial differential equation 
• associated with a Riemann-Liouville derivative of fractional order as the time-derivative and Riesz-Feller 

fractional derivative as the space derivative. The method followed in deriving the solution is that of joint 
Laplace and Fourier transforms. The solution is derived in a closed and computable form in terms of the 
H- function. It provides an elegant extension of the results given earlier by Debnath [4], Chen et al. [2], 
Haubold et al. [17], Mainardi et al. [27,28 ], Saxena et al. [39], and Pagnini et al. [35]. The results obtained 
are presented in the form of four theorems. Some results associated with fractional Schrodinger equation 
and fractional diffusion-wave equation are also derived as special cases of the findings. 
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ON ' 1. Introduction 

o . 

General models for reaction-diffusion systems are discussed by Henry and Wearne [18, 19], Henry et 
al. [20], Haubold et al. [16,17], Saxena et al. [38, 39, 40], Mainardi et al. [27, 28] and others. Stability 
in reaction-diffusion systems and nonlinear oscillation have been discussed by Gafiychuk et al. [12,13]. 
Recently, Engler [7] obtained the speed of spread for fractional reaction-diffusion. Distributed order sub- 
diffusion is discussed by Naber [32]. In a recent paper, Chen et al. [2] have derived the fundamental and 
numerical solution of a reaction diffusion equation associated with the Riesz fractional derivative as the 
space derivative. Reaction-diffusion models associated with Riemann-Liouville fractional derivative as the 
time fractional derivative and Riesz-Feller derivative as the space derivative are recently discussed by Haubold 
et al. [17]. Such equations in case of Caputo fractional derivative are also solved by Saxena et al. [39]. In 
connection with the evolution equations for the probabilistic generalization of the Voigt profile function, it 
is shown by Pagnini and Mainardi [35] that the solution of the following integro-differential equation 

ON 

— = D« 1 N(x,t)+ D« 2 N(x,t), N{x,0) = S(x) (1.1) 

in terms of its Fourier transform, where qDq 1 and qDq 2 are the Riesz fractional derivatives of orders ai and 
«2 respectively, and S(x) is the Dirac-delta function, which is given in [35, p. 1593]. Consider the following 
Fourier transform, with the Fourier parameter k: 

N*(k,r) = exp{-T(|A:| ai + \k\ a *)}. (1.2) 
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This has motivated the authors to investigate the solutions of fractional partial differential equations (3.1) 
and (3.14). The technique used in deriving the results is that of joint Laplace and Fourier transforms. The 
results are obtained in a closed and computable form. Due to the general character of the derived results, 
many known results given earlier by Chen et al. [2], Haubold et al. [17], Saxena et al. [39], Pagnini et al. 
[35], and others readily follow as special cases of our derived results. Solutions of certain extended space-time 
fractional diffusion wave equation and generalized Schrodinger equations are also deduced from our findings. 

The Schrodinger equation is a fundamental equation of quantum mechanics. Feynman and Hibbs [11] re- 
constructed the Schrodinger equation by making use of the path integral approach by considering a Gaussian 
probability distribution. This approach is further extended by Laskin [23-25] in formulating the fractional 
Schrodinger equation by generalizing the Feynman path integrals from Brownian-likc to Levy-like quantum 
mechanical paths. In a similar manner, one obtains a time fractional equation if non-Markovian evolution 
is considered. In a recent paper, Naber [32] discussed certain properties of the time fractional Schrodinger 
equation by writing the Schrodinger equation in terms of fractional derivatives as dimensionless variables. 
The time fractional Schrodinger equations are discussed by Debnath [4,5], Bhatti and Debnath [1], Saxena 
et al. [38], Tofight [40], and others. 

2. Mathematical Prerequisites 

The Riemann-Liouville fractional derivative of order a > is defined as (Samko et al. [37, p. 37]; also 
see, Kilbas et al. [21]) 

i c m ft 

D?f(x, t) = - — / (t - r)"-"- 1 /^, t)dt, n = [a] + 1, n e N, t > 0, (2.1) 

where [a] means the integer part of the number a. The Laplace transform of the Riemann-Liouville fractional 
derivative is given by Oldham and Spanier [34, eq.(3.1.3)]; (also see, Kilbas et al. [21]): 

n 

L{ D?N(x,t);s} = s a N{x,s) -^s r - 1 aDf- r N(x,t)\t=o, n - 1< a < n. (2.2) 

r=l 

This derivative is useful in deriving the solutions of integral equations of fractional order governing certain 
physical problems of anomalous reaction and anomalous diffusion. In this connection, one can refer to the 
monograph by Podlubny [36], Samko et al. [37], Oldham and Spanier [34], Miller and Ross [31], Kilbas et 
al. [21], Mainardi [26], Diethelm [6], and recent papers on the subject [27,28,30,33,40,41]. Following Feller 
[9,10], it is conventional to define the Riesz-Feller space-fractional derivative of order a and skewness 9 in 
terms of its Fourier transform as 



F{ x D%f(xy,k} = -* e a (k)f*(k), (2.3) 

where 

9ir 

\f> e a (k) = \k\ a e*v[i{signk)—), < a < 2, \6\ < min{a, 2 - a}. (2.4) 
Further when 9 — 0, we have a symmetric operator with respect to x that can be interpreted as 

x DS = -{-^y. (2.5) 

This can be formally deduced by writing —(k) a = — (k 2 )% . For < a < 2 and \9\ < min{a, 2 — a}, the 
Riesz-Feller derivative can be shown to possess the following integral representation in x domain: 

.WW = ^{ Mia + 9 W] f ^±g^ d C 

+ sin[(a - jf° f{x f{x) d(} (2.6) 
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For 9 — 0, the Riesz-Feller fractional derivative becomes the Riesz fractional derivative of order a for 
1 < a < 2, defined by analytic continuation in the whole range < a < 2, a / 1 (see, Gorenflo and 
Mainardi [14]) as 



X D% = -A[J+ a - IZ a ], 



where 



A = 



1 H 2 

* . j-a _ u j2-q 

2 cos(a7r/2) ' + _ Ax 2 + 



The Weyl fractional integral operators are defined in the monograph by Samko et al. [37] as 



(2.7) 
(2.8) 



and 



1 f°° 

(I P + N)(x) = —jj x - C)^ 1 iV(C)dC, > o 
1 f°° 

(I P -N)(x) = YiJ ) j x (C - xf- l N{QdC, (3>0. 



Note 1. We note that x Dq is a pseudo differential operator. In particular we have 



,Dl = , but x Dl * A. 



(2.10) 



For 8 = 0, we have 



F{ x D%f(xy,k} = -\k\ a f*(k). 



(2-11) 



The H-function is defined by means of a Mellin-Barnes type integral in the following manner [29, p. 2]: 



H p,q ( Z ) - H p 



m,n 

q 



Trm,n 
P,Q 



\(a p ,A p ) 
\{b q ,B q ) 

\(ai,A 1 ),...,(a p ,A p ) 
l(6i,Bi),...,(6 9 ,B 9 ) 



where i — (—1) 2 , 



6(0 = 



{Il?=i £fe + ^OHTI^i r(i - oj - A,-Q> 
fflUi r(i - b 3 BjOHUU+i r (°i + ^0} : 



(2.12) 



(2.13) 



and an empty product is always interpreted as unity ; m, n, p, q E No with < n < p, 1 < m < q, 
Aj,R+, dj,bj € R or C ', i = l...,p; j = 1, q) such that 



-Ai (6j + fc) ^ Bj(a, - s - 1), fc,s e iV i = 1, ...,n; j = 1, ...,m, 



(2.14) 



and these poles are separated, where we employ the usual notations: N a = (0,1,2,...); R = (—00,00), 
and C being the complex number field. A comprehensive account of the H-function is available from the 
monographs Mathai et al. [29] and Kilbas et al. [22]. We also need the following result in the analysis that 
follows; Haubold et al. [16] has shown that 



F- 1 [E M {-at^ e a {k);x] = —H, 



2,1 
3,3 



(2.15) 



where Sft(a) > 0, &(/?) > 0, 5R( 7 ) > 0. 

3. Solution of Unified Fractional Partial Differential Equations 
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In this section, we will investigate the solution of fractional partial differential equations, which may be 
regarded as an extension of one-dimensional fractional reaction-diffusion equation, one-dimensional space- 
time fractional diffusion-wave equation and one-dimensional fractional Schrodinger equation. The results 
are presented in the form of the following four theorems. 

Theorem 1. Consider the following one- dimensional non-homogeneous unified fractional differential equa- 
tion: 

a D?N(x, t) = A x DgN(x, t) +f i x DlU(x, t), (3.1) 
where t > 0, x G R; a, 6, (3, 7 and <p are real parameters with the constraints 

< [3 < 2, < 7 < 2, \6\< min(/3, 2 - B), < min( 7 , 2 - 7), < a < 1, (3.2) 
with the initial conditions 

o£> t Q_1 ^V(x,0) =f(x), forx e R, lim N(x, t) = 0, t > 0. (3.3) 

X -^±00 

Here £)" _1 Ar(x, 0) means the Riemann-Liouville fractional partial derivative of N(x,t) with respect to t of 
order a — 1 evaluated at t = 0; x Dg and are i/ie Riesz-Feller space-fractional derivatives respectively of 
orders [3 and 7 wif/i asymmetries 9 and <j>. n D" is the Riemann-Liouville time-fractional derivative of order 
a; A and /i are arbitrary constants, and f(x) and U = U(x,t) are the given functions. Then for the solution 
of (3.1 ), subject to the above conditions, there holds the formula 



G(x-U)f(0^ 

-00 



(i / (t-rr^i G 1 {x-T,t-T)U{^T)A^]dT, (3.4) 



where the Green functions G{x,t) and Gi(x,t) are given by 

j-a-l r°° 



G(x,t) = / exp(-ikx)E aa (-X^ e (k)t a )dk, (3.5) 

2?r J-oc 

(3.6) 



m 33 



(A t«)> 



(l,^),(a,f),(l,p)- 



(1,1), (1 ,£),(!,*>) 



wii/i /3 > 0, p = and 

1 Z" 00 

G 1 (x,t) = — exp(-ikx)9+(k)E ata (-\ ^(k)t a )dk (3.7) 
2 71 " J -00 

and H 3 ' 3 (z) is the H-function defined by (2.12) and E a> p(z) is the Mittag-Leffler function, defined by [8, 
Section 18.1] 

E a . ( 3(z) = V — =r, a, 6 e C, »(a) > 0, 5R(/3) > 0. (3.8) 

^ r(an + /3) 

In deriving the value of the integral in (3.5), the formula (2.15) has been used. 
Proof: In order to derive the solution of (3.1), we introduce the joint Laplace-Fourier transform in the form 



*■ 00 />OC 

-st+ikx 



a 



N*(k,s) = / c- st+lKX N(x,t)dx dt, (3.9) 



where 3fJ(s) > 0, k > 0. If we apply the Laplace transform with respect to the time variable t, Fourier 
transform with respect to space variable x and use the initial conditions (3.2), (3.3) and the formula (2.2), 
then the given equation transforms into the form 



s a N*{k, s) - f*(k) = -X V e a {k)N*{k, s)-n ^{k)U*{k, s), 

where according to the conventions followed, the symbol N(k,s) will stand for the Laplace transform with 
respect to time variable t and * represents the Fourier transform with respect to space variable x. Solving 
for N*(k,s), it yields 



N*(k,s) = 



f*(k) fx^{k)U*{k,s) 



(3.10) 



s« + A *«(£;) s a + \^ f} {k) 

To invert (3.10), it is convenient to first invert the Laplace transform and then the Fourier transform. Thus 
to invert the Laplace transform we use the formula given in [39] 



a + s a 



(3.11) 



where 5ft(s) > 0, $l(a) > 0, 3ft(a — (3) > —1 and the convolution theorem of the Laplace transform to obtain 



N*(k,t) = r{k)t a - i E a>a {-\ n^(k) 

-t" *t(k) f (t - ^^(-a^Mti - T) a )U*{k,T)&T. 



(3.12) 
Jo 

Now the application of the inverse Fourier transform gives the exact solution in the following form: 

N(k,t) =t a - 1 F- 1 [f*(k)E a , a (~\ ^(k)t a )} 

-»[ {t-Tr^F-^^Ea^-X^kXt-T^U^k^dT. (3.13) 
Jo 

Finally, the application of the convolution theorem of the Fourier transform yields the desired solution (3.4). 
Following a similar procedure, it is not difficult to establish the following theorems: 

Theorem 2. Consider the same equation in (3.1) with the same condition on the parameters except that 
< a < 2, instead of < a < 1. In addition to the initial conditions in (3.3), assume that 



D?- 2 N(x,0)=g(x). 
Then the solution of (3.1) under the conditions (3.2), (3.3) and (3.14) is given by 



(3.14) 



G(x - C t)/(£)d£ + / G 2 {x-Z, t)g{Od$ 

- oo J — OO 

-H {t-r) a - l [ G 1 (x-r 1 t- T )U(tr)dS]AT, 

JO J -oo 



(3.15) 



where the Green functions G(x,t) and G\{x,t) are defined in (3.6) and (3.7) respectively and G2(x,t) is 
given by 



j-ot — 2 roo 

G 2 (x,t) = — — / exp(-ikx)E a ^ 1 (-\ Vt(k)t a )dk 

271" J-oo 

a-2 . T Ul (l.?).(«-l.f).(l.P)l 



t 



f)\x\ H ™ 
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(i,i),(i,£),(M 



(3 > 



(3.16) 



with p — ^2]T and H 3 ' 3 (z) is the H-function defined by (2.17) and E a< p(z) is the Mittag-Leffler function, 
defined in (3.8). 

Now, we set U(x,t) — N(x,t), where N(x,t) is the unknown function to arrive at the next result. 

Theorem 3. Consider the same equation in (3.1) with U(x,t) replaced by N(x,t). Assume that the condi- 
tions (3.2) and (3.3) hold. Then solution of (3.1) with U(x,t) replaced by N(x,t) is given by 

/oo 
G 3 (x-Z,t)f(Odt, (3.17) 
-oo 

where the Green function Gs(x,t) is given by 

j.a—1 r oo 

G 3 (x,t) = — — / exp(-ikx)E a , a [-(\ *«(fc) + M ^(k))t a }dt. (3.18) 



Similarly, we find that the following theorem holds true: 

Theorem 4. Consider the same equation in Theorem 3 with < a < 2 and under the conditions (3.2), 
(3.3), and (3.14). Then the solution of (3.1) with U(x,t) replaced by N(x,t) with < a < 2 is given by 

/oo 
G 3 (x-C,t)f(OdL (3.19) 
-co 

where the Green function Gs(x,t) is defined in (3.18) and the other Green function G^{x,t) is given by 

+a-2 roc 

Gi{x,t) = — — / exp(-ifca;)£ Q , Q _i[-(A ^(k)+p^(k))t a ]dk. (3.20) 

2 71 " J-oc 

Note 2. It is interesting to observe that for g(x) — 0, Theorems 2 and 4 yield similar types of results as 
Theorems 1 and 3 respectively. 

4. Selected Special Cases 

If we set 9 — cf) = then by virtue of the identity (2.11), Riesz-Feller derivative reduces to Riesz 
derivative and consequently Theorem 1 yields the following results: 

Corollary 1. Consider the one-dimensional non-homogeneous unified fractional differential equation (3.1) 
for 6 = <f> = 0: 

D<?N{x, t) = A x D$N(x, t) + fx x DlU{x, t), (4.1) 

where the conditions (3.2) and (3.3) hold under 9 = <fi = 0, t > 0, x G R, a, (3, 7 are real. Then for the 
solution of (4-1) there holds the formula 



N(x,t) = J 



-M/ {t-rT^F-^ Ge(x-T,t-T)U(Z,T)dZ]dT, 

JO J -00 



(4.2) 



where the Green functions G§{x,t) and Ge{x,t) are given by 

j-a-l roo 



j.a—1 roo 

G 5 (x,t) = - / exp(-ikx)E aa (-X\kft a )dk (4.3) 

(l,i),(a,f),(l,i)- 
(l,l),(l,£),(l,i) . 
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(3 > (4.4) 



i r°° 

Ge(x,t) = — cxp(-ifc.x)|fc| 7 S Q i(-A|fc| /3 t")dfc. (4.5) 

2lT J-oo 



For fi = <f> = 0, Theorem 1 reduces to the following result given by Haubold et al. [10]. 
Corollary 2. Consider the following one- dimensional fractional reaction- diffusion model of (3.1) with fi = 

o D?N(x,t) = \ x D e N(x,t), (4.6) 

where t > 0, x e R, a, 9, (3 real, and conditions (3.2) and (3.3) hold under p, = <fi = 0. Then for the 
solution of (4-6) there holds the formula 

/OO 
G(x - Z,t)f(Z)dt, (4.7) 
-oo 

where the Green function G(x,t) is defined in (3.6). 

If we further set f(x) = S(x), where 5(x) is the Dirac-delta function, we obtain the fundamental solution 
of the space-time fractional diffusion equation given by Haubold et al. [17]. 

By setting 9 = <f> = and using (2.11), the following Corollaries 3-5 can be easily deduced from Theorems 

2-4. 

Corollary 3. Consider the following one- dimensional non-homogeneous unified fractional differential equa- 
tion: 

D?N(x, t) = A x D$N(x, t) + n x D^ U(x, t), (4.8) 
where (3.2), (3.3), (3.14) hold under 9 = <f> = 0. Then for the solution of (4-8), there holds the formula 

/>oc />oo 
G 5 (x - £, t)f{Z)d$ + / G 7 (x - C t)g{S)d$, 
-oo J — oo 



-H {t-T) a ^F-\ G 6 (x-T,t-T)U(Z,T)d£]dT, (4.9) 
JO J-oo 

where the Green functions Gs{x,t) and Ga(x,t) are defined in (4-3) and (4-5) respectively and Gr(x,t) is 
given by 

j-a-2 /■oo 

G 7 (x, t) = / cxp(-ikx)E a a (~\\kft a )dk 

2tt J^oo 

> 0, (4.10) 



j.a-2 

l rr2,l 

f}\x\ " 3 < 3 



(At a )1 



JW 1 ' 2 ) 



where p = . 



Corollary 4. Consider the following one- dimensional non-homogeneous unified fractional differential equa- 
tion where a, f3 and 7 are real parameters with the constraints 

oD?N(x, t) = A x D$N(x, t)+fi x DlN{x, t), (4.11) 
where (3.2), (3.3) hold for9 = cj> = 0. Then 
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/oo 
G 8 (x - £,t)f(Qd$, (4.12) 
-oo 

where the Green function G$(x,t) is given by 

J-Q-l /■OO 

G s (x,t) = / exp(-ikx)E aa [-(\\kf +n\k\~<)t a }dk. (4.13) 

2tt ./-oo 



Corollary 5. Consider the following one- dimensional non-homogeneous unified fractional differential equa- 
tion 

D?N(x, t) = A aDgiVfo t) + fi x DlN{x, t), (4.14) 
where (3.2), (3.3), (3.14) hold. Then 

/OO p OO 

Ga(x-S,t)f(t)d$+ t G 9 (x - Z,t)g(t)d$, (4.15) 
-oo J —oo 

where the Green functions Gs{x,t) is defined in (4-13) and the other Green function Gg(x,t) is given by 

j-a-2 roo 

G 9 (x,t) = —— cxp(-ikx)E a . a ^[-(X\kf + ^\ky)t a ]dk. (4.16) 

27r J-oo 

Note 3. It is observed that for g(x) = 0, Theorems 2 and 3 yield similar types of results as Theorems 1 
and 3, respectively. 

5. Further Special Cases 

In this section we discuss some consequences of the main results relating to the fractional Schrodinger 
equation. In the following corollaries 6 and 7, h is the Planck constant. 

If we set 7 — <f> = and A = then Theorems 1 and 2 give the solutions of non- homogeneous fractional 
Schrodinger equations as shown below. 

Corollary 6. Consider the following one- dimensional non-homogeneous unified fractional Schrodinger 
equation of mass m: 

a D?N(x, t) = (^-) x D%N(x, t)+n U(x, t), (5.1) 
where t > 0, x e R, a, [3 real, with the constraints 

< /3 < 2, < a < 1, (5.2) 

with the initial conditions (3.3), and N(x,t) as the wave function. 

Then for the solution of (5.1 ), subject to the above constraints, there holds the formula 

/oo 
G(x-U)f(m 
-OO 

/OO pOC 
(t-O^i Gio(x-S,t-T)U(Z,T)d£]dT, (5.3) 
-oo J — oo 

where the Green function G(x ) t) is defined in (3.6) and G\o{x,i) is given by 

/CO 
exp(-ikx)E ata [-(\ 1> 9 p (k))t a }dk 
-co 



1 u 2 - 1 

3\x\ s ' 3 



(at a )> 



(l,£),(a,f),(l,P) 
(l,l),(l,£),(l,p) 



(5.4) 



where (3 > 0, a = ^, P=^- 

Corollary 7. Consider the following one- dimensional non-homogeneous unified fractional Schrddinger 
equation of a particle of mass m: 

D ?N(x, t) = (^-) x D%N(x, t)+n U(x, t), (5.5) 
where t > 0, x e J?, a, 0, /3 reaZ 7 wit/i i/ie constraints 

< /3 < 2, |0| < min(/3, 2-/3), 1 < a < 2, (5.6) 

and with the initial conditions (3.3). Then for the solution of (5.5), subject to the above constraints, there 
holds the formula 

/oo /.oo 
G(x - £, t)f(Od£ + tG 2 (x-£, t)g(Od£ 
-oo J — oo 

ft pOO 

+ V / (t-rr^i Gio(x-T,t-T)U(Z,T)d£]dT, (5.7) 

JO J-oo 

where the Green functions G(x,t),G2(x,t) and Gio(x,t) are respectively given by (3.6), (3.16), and (5.4), 
respectively. 

On the other hand, if we set 7 = = 0, Theorems 1 and 2 provide the solutions of non-homogeneous 
fractional generalized diffusion-wave equations as given below: 

Corollary 8. Consider the following one- dimensional non-homogeneous unified fractional diffusion-wave 
equation: 

oD?N(x, t) = c 2 x D t3 g N{x, t)+fi U(x, t), (5.8) 
where t > 0,x G R, a, j3 real, with the constraints 

0</3<2, < a < 1, (5.9) 

and with the initial conditions (3.3) and (5.2), where fi and c are arbitrary constants. Then for the solution 
of (5.8), subject to the above conditions, there holds the formula 

/oo 
G(s-£,t)/(Od£ 
-OO 

/OO /* 00 

G w {x-tt-T)V%T)d£\&T, (5.10) 
-00 J —oo 

where the Green function G(x,t) and G\o{x,t) are respectively given in (3.6) and (5.4). 

Corollary 9. Consider the following one- dimensional non-homogeneous unified fractional diffusion-wave 
equation: 

oD?N(x, t) = c 2 x D p g N(x, t)+n U(x, t), (5.11) 

where t > 0, x <E R, a, 9, (3 real, with the constraints < (3 < 2, \9\ < min(/3, 2 — !3), 1 < a < 2, and with 
the initial conditions (3.3) and (3.14). Then for the solution of (5.11), subject to the above conditions, there 
holds the formula 

9 



/OO POO 
G(x - £, i)/(Od£ + / tG 2 (x-£, t)g(Od£ 
-oo J — oo 

ft pOO 

+ H / (i-r)"- 1 !/ Gio(x-T,t-T)^,T)d^]dT, (5.12) 

JO ./-oo 

where the Green functions G(x,t), G2(x,t) and G\o(x,t) are, respectively, given by (3.6), (3.16), and (5.4). 

Now, setting fi = 1 and 9 = (3 = 0, we find that the Corollaries 8 and 9 give rise to the expressions for the 
solution of non-homogeneous one-dimensional fractional generalized wave equations as given by Debnath [4, 
p.141]. 

Finally, for a = 1 and 6 = <f> = 0, Theorem 3 provides a solution of generalized diffusion equation with two 
space-fractional derivatives, which were recently studied by Pagnini and Mainardi [35]. 
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